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ABSTRACT 


Wt»~Gtudy  ^transonic  flows  along  a  nozzle^ based  on  a  one-dimensional 
model.  It  is  shown  that  flows  along  the  expanding  portion  of  the  nozzle  are 
stable.  On  the  other  hand,  flows  with  standing  shock  waves  along  a 
contracting  duct  are  dynamically  unstable.  This  was  conjectured  by  the  author 
based  on  the  study  of  noninteracting  wave  patterns.  The  author  had  shown 
earlier  that  supersonic  and  subsonic  flows  along  a  duct  with  various  cross 
sections  are  stable.  Basic  to  our  analysis  are  estimates  showing  that  shock 

waves  tend  to  decelerate  along  an  expanding  duct  and  accelerate  along  a 

i 
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contracting  duct, 


AMS(MOS)  Subject  Classifications:  76H05,  76E30,  35L65,  35L67 

Key  Words:  Transonic  flows,  stability  and  instability,  flows  along  a  nozzle, 

asymptotic  behavior. 

Work  Unit  Number  1  -  Applied  Analysis 

c 

“7' 


Department  of  Mathematics,  University  of  Maryland,  College  Park,  MD  20742. 


Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 .  This 
material  is  based  upon  work  supported  by  the  National  Science  Foundation  under 
Grant  No.  MCS  7802202  and  by  the  Sloan  Foundation. 


NTIS  GRA&I  ~~gj 
DTIC  TAB 
Unannounced 
Justification.. 


By - 


STABILITY  OF  TRANSONIC  FLOWS  ALONG  A  NOZZLE 


Distribution/ 


Availability  Codes 

Avail  and/or 
Dist  Special 


Tai-ping  Liu 

$1.  Introduction 

It  is  well-knovm  that  gas  flows  are  often  highly  unstable.  This  is  clear  from 
experimental  studies  and  causes  great  difficulties  in  numerical  calculations.  However,  the 
emergence  of  shock  waves  in  the  flow  as  a  consequence  of  the  nonlinearity  of  the  gas 
dynamics  equations  makes  it  even  more  difficult  for  analytical  studies.  The  present  paper 
is  a  step  in  understanding  such  phenomena.  We  show  that  for  a  1-dimensional  model  flows 
along  an  expanding  duct  are  always  asymptotically  stable)  while  flows  with  a  standing  shock 
wave  along  a  contracting  duct  are  dynamically  unstable.  Since  the  initial  value  problem 
for  the  gas  dynamics  equations  is  expected  to  be  well-posed,  by  instability  we  mean'  the 
asymptotic  state  of  a  flow  does  not  depend  smoothly  on  the  initial  state. 

The  equations  which  model  gas  flow  in  a  variable  area  duct  are  of  the  form  ((1)): 


(1.1) 
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where  a(x)  is  the  cross  section  of  the  duct,  p,  u,  p  and  E  are,  respectively,  the 
density,  velocity,  pressure  and  the  total  energy  of  the  gas,  and  x  e  R,  t  >  0.  For  a 
uniform  duct,  a'(x)  =  0,  the  equations  become 
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In  this  case  it  was  shown  by  Liu  [6]  that  the  asymptotic  state  of  a  solution  of  the  initial 
value  problem  consists  of  the  elementary  waves  determined  by  the  values  of  the  initial  data 
at  x  «  ±°».  In  particular,  the  flows  are  always  asymptotically  stable.  The  reason  for 
such  a  strong  stability  phenomenum  is  that  in  spite  of  very  complicated  nonlinear  wave 
interactions,  waves  also  tend  to  combine,  cancel  and  decouple  as  a  result  of  the  dependence 
of  the  characteristic  speeds  on  the  state  and  the  hyperbolicity  of  the  system. 

Consequently,  a  general  solution  tends  to  a  wave  pattern  which  is  noninteracting.  Such  a 
wave  pattern  consists  of  shock  waves,  rarefaction  waves  and  a  traveling  wave,  and  can  be 
found  by  solving  the  Rlemann  problem  which  is  an  Initial  value  problem  with  the  Initial 
data  containing  only  two  constant  states: 

for  x  <  0  , 
u>r  for  x  >  0  , 

(i)  2  ( p ,  pu ,  pE ) 


where  w.  and  in  are  constant  states, 
i  r 

For  variable  cross  section  duct,  it  has  been  shown  by  Liu  [7]  that  supersonic  and 
subsonic  flows  are  also  asymptotically  stable  provided  that  the  duct  becomes  uniform  as 
x  ♦  ±*».  In  this  case,  the  asymptotic  state  consists  of  the  aforementioned  elementary 
waves  and  a  standing  wave  which  is  a  solution  of 
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When  the  flow  ts  not  transonic/  shock  waves  and  rarefaction  waves  travel  with  nonzero  speed 
and  eventually  decouple  from  the  standing  wave.  For  transonic  flows,  these  two  kinds  of 
nonlinear  modes  intertwine  and  more  interesting  nonlinear  wave  phenomena  occur.  Based  on 
the  above  studies,  one  expects  that  a  general  solution  also  tends  to  an  asymptotic  state 
which  does  not  create  any  nonlinear  interaction.  Such  noninteracting  wave  patterns  were 
studied  in  Liu  (9].  It  was  found  that  for  an  expanding  duct  the  wave  pattern  is  uniquely 
determined  by  its  value  at  x  «  ±«  and  thus  one  should  expect  the  general  solution  to  be 
asymptotically  stable.  On  the  other  hand,  along  a  contracting  duct  there  may  exist  three 
noninteracting  wave  patterns  with  same  given  end  states  at  x  »  ±°°.  one  wave  pattern  is 
smoothly  deformed  to  a  supersonic  flow  and  the  other  to  a  subsonic  flow.  Thus  these  two 
should  be  stable  by  [8].  The  third  one,  which  contains  a  standing  shock  wave,  is  therefore 
expected  to  be  unstable  by  the  criterion  of  the  exchange  of  stability.  It  is  the  purpose 
of  this  present  paper  to  verify  the  above  conjecture  on  the  stability  of  wave  patterns  in 
the  dynamical  sense.  In  particular,  we  show  that  when  a  standing  shock  wave  in  a 
contracting  duct  is  perturbed,  it  accelerates  itself  and  eventually  tends  to  a  stable  wave 
pattern  with  no  standing  shock  wave  which  differs  greatly  from  the  original  unperturbed 
state.  On  the  other  hand,  when  a  standing  shock  wave  in  an  expanding  duct  is  perturbed,  it 
decelerates  and  tends  to  a  nearby  asymptotic  state.  Actually,  as  waves  move  along  the  duct 
complicated  wave  interactions  begin  to  occur,  nevertheless,  the  above  statements  remain 
valid  qualitatively. 

Our  basic  estimate  is  on  the  change  of  speed  of  a  shock  wave  as  it  moves  along  a 
duct.  This  and  other  estimates  on  wave  interactions  are  established  in  Section  3.  Our 
analysis  is  based  on  the  random  choice  method  of  Liu  [8] ,  which  generalizes  the  Glimm 
scheme,  [3] ,  for  conservation  laws.  The  idea  is  to  decompose  the  solution  into  nonlinear 
inodes  for  conservation  laws  and  standing  waves.  In  Section  4  we  describe  a  variation  of 
the  simplified  method  of  Liu  [9J.  Finally,  in  Sections  S  and  6,  respectively,  we  study 
flows  along  a  contracting  duct  and  along  an  expanding  duct. 


$2.  Preliminaries 

The  gas  dynamics  equations  (1.2)  are  supplemented  by  the  constitutive  relation 

p  -  p(v,s) 

where  s  is  the  entropy  of  the  gas,  v  -  1/p  the  specific  volume.  These  physical 
quantities  satisfy  the  thermodynamics  relation  de  ■  Tds  -  pdv,  T  the  temperature  of  the 
gas.  He  assume  that 


3p/ 3v  <  0  and  32p/  2  >  0 

dV 


The  characteristic  speeds  X^»  i  ■  1,2,3,  and  corresponding  characteristic  right  vectors 
r,  are  (2J : 


V,  v, 

X1  »  u  -  (-  3p/3v)  2v,  X2  “  u,  Xj  »  u  +  (-  3p/3v)'2' 


Y,  ”  ( 1 ,  X , , E  +  pv  -  uv{-  3p/3v)/'2  ) )  , 


Y3  ■  ( 1 ,  Xj,E  +  pv  +  uv(-  3p/3v)'2))  , 


and  y2  is  such  that  y2  •  Vu  ■  y2  •  Vp  »  0.  Hereafter  the  coordinate  system  for  the 
states  is  ui  »  (p,pu,pE).  The  first  inequality  in  (2.1)  implies  that  X^,  1  «  1,2,3,  are 
real  and  distinct  and  so  systems  (1.1)  and  (1.2)  are  strictly  hyperbolic.  The  second 
inequality  in  (2.1)  implies  that 


Y,  •  7X,  -  -  V2  32p/3y2  ♦  (-  3p/3v)"  1/2v3/2  <  0  , 
Y3  •  ’Xj  -  V2  32p/3v2  •  {-  3p/3v)' 1/zv3/2  >  0  . 
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Direct  calculations  also  yield: 


(2.5)  Yl  .  Vu  =  -  (-  ap/avj^v2  <  0.  f3  •  Vu  -  (-  3p/3v)1/2v2  >  0  , 

2 

(2.6)  y  *  Vp  “  y3  •  7p  “  -v  3p/3v  >0  . 

We  now  describe  the  elementary  waves  for  Euler  equations  (1.2).  An  i-rarefactlon  wave 
(  <dq  ,  (D.J ) ,  i  *  1,2,3,  is  a  smooth  solution  with  speed  X^  which  depends  on  one  parameter 
and  connects  uiQ  on  the  left  and  on  the  right.  It  is  a  consequence  of  (1.2)  that 

is  on  the  i-rarefaction  curve  R^(Ug),  the  integral  curve  of  r^  through  ,  and 
Xi  ( U) )  >  X.tV: 


£  R*<V 


{<u  e  Rt  ( ( to)  >  (  ojq  ) }  . 


There  are  two  kinds  of  discontinuity  waves.  For  the  2-characteristic  field,  ( ,  ti)^ )  is  a 
discontinuity  wave  with  speed  0(^,0^)  provided 


£  T2  ( S  !  =  u(u0),  p(  u>)  *  pfiiip)} 


and  is  a  constant  discontinuity: 


<J<  ,  u1 ) 


V  V 


V  V 


For  the  first  and  third  characteristic  fields,  a  discontinuity  wave  (u)g,w^)  with  speed 
oIwq.w.j)  satisfies  the  following  jump  ( Rankie-Hugoniot)  condition: 


and  is  a  stock  wave  satisfying  the  following  stability  property,  [1): 


( oi1 )  <  o(  ii)0 ,  { u»Q ) 


If  we  set 


T^uip)  =  Si<“o)  u  1  ”  1,3  ' 


then  any  state  u>  on  T^(uQ}  can  be  related  to  on  the  left  by  an  i-elententary  wave 


(u,p)-plane: 


(2.7) 
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The  above  holds  in  a  small  neihborhood  of  w  under  the  assumption  (2.1).  It  holds 
globally  under  a  stronger  assumption,  (5): 


(2.8) 


p  -  p(v,e),  3p/3v  <  0,  3p/pe  >  0  . 


It  follows  from  the  thermodynamics  relation  de  =*  Tds  -  pdv  that  3p/3e  -  T  3p/3s  and  so 
the  second  inequality  in  (2.8)  is  equivalent  to,  [11]: 

p  »  p(v,s) ,  3p/3s  >  0  . 


To  solve  the  Riemann  problem  (1.2)  and  (1.3)  we  find  states  u>  and  u>  such  that 

m  n 

u  e  T^(o),),  (ii  e  T_  ( (i)  )  and  e  T,(u  )  so  that  the  solution  consists  of  1-wave 
m  i  Jt  n  z  m  n  3  r 

(u.,(u  ),  2-wave  (a  ,u>  >  and  3-wave  ( a>  ,01  ).  The  states  to  and  to  are  determined 
im  mn  nr  mn 

uniquely  by  to^  and  due  to  (2,7).  We  choose  any  nonsingular  parameter  t^, 

i  «  1,2,3,  along  T^  curves  to  measure  the  strength  of  waves: 


(vVi  =  W  •  W 


‘VVa  S  W  '  T2 ' un ’  ' 


S'Va  *  T3 ‘ “n ^  _  W  ‘ 


We  choose  so  that  shock  waves  have  negative  strength  and  rarefaction  waves  have 

positive  strength. 


a 

*:> 


I 


i 

» 

«*■ 

i 
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$3.  Basic  estimates 

The  steady  state  equations  (1.4)  can  also  be  written  as 


.  r  («)  r  (oil 

^»V2c(x)pu(— +  — )  . 


where  c(x)  =  -a'(x)/a(x)  which  is  assumed  to  be  small  in  absolute  value.  In  this  section 
we  will  fix  two  space  positions  x  *  x_,  x  •  x+(  x_  <  x+,  and  call  (<i>q,u.|)  a  standing 
wave  if  there  exists  a  solution  u(x)  of  (1.4),  or  equivalently  (3.1),  such  that 
<u(x  )  *  (dq  and  <*i(x+)  ■  We  will  consider  transonic  flows  with  positive  speed,  that 

is,  (  X 1  f  to)  I  is  small.  We  first  study  the  evolution  of  subsonic  and  supersonic  waves  as 
they  propagate  along  the  duct. 

Lemma  3.1.  Suppose  that  all  the  states  are  not  sonic,  that  is,  there  exists  \t  >  0  such 
that  either  X^(w)  >  X.  for  all  u  under  consideration  or  X^(id)  <  -X,  for  all  m  under 
consideration.  Let  (oi^.u^)  and  (u^,^)  be  two  standing  waves,  then  the  strengthes  of 
elementary  waves  in  the  solutions  of  the  Riemann  problems  (u^,^)  and  (01^,111^.)  for  (1.2) 
are  related  as  follows: 

‘VVi  =  *  V  Vi  +  0(1)<jc(X«)‘2,  i  -  1,2,3  , 

a  S  , 

c  5  max  |c(x)|  , 

x<x<x+ 


where  0(1)  is  a  bound  depending  only  on  the  system  (1.2). 

Proof:  Since  the  waves  T^,  i  =  1,2,3,  are  smooth,  the  solution  of  the  Riemann  problem 
depends  continuously  on  the  end  states.  Thus  the  lemma  follows  from  the  elementary  theory 
of  ordinary  differential  equations.  We  include  (X,)  on  the  right  hand  side  of  the 
estimate  just  to  remind  us  that  X#  may  be  a  small  number. 

(J.E.D. 
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The  following  crucial  estimate  is  on  the  evolution  of  a  1-shock  wave  which  is  a 
perturbation  of  a  standing  1-shock  wave.  It  shows,  in  particular,  that  the  shock  wave 

accelerates  as  it  propagates  along  a  contracting  duct  and  decelerates  as  it  propagates 

along  an  expanding  duct. 

Lemma  3.2.  Suppose  that  (10^,10^)  is  a  1-shock  wave  which  is  a  perturbation  of  a  standing 

shock  wave,  that  is,  X^(io^)  >  0  >  X^fuO  and  its  speed  <j  5  0(10^, u^,)  is  small  compared 

to  its  strength  a  =  I  <  to^,  io^ )  I .  Suppose  also  that  the  duct  is  not  uniform  for 

x  c  (x  ,x+),  in  other  words,  either  c/2  <  c(x)  <  c,  x_  <  x  <  x+  or  c/2  <  -c(x)  <  c, 

x  <  x  <  x  ,  for  some  small  positive  constant  c.  Let  (io.,io.)  and  (to  ,io  )  be  two 

+  l  l  r  r 

standing  waves,  then  there  exist  to.  on  S^(io^)  and  a  positive  constant  d  depending 
only  on  (1.2)  such  that 
(i)  II oo.  -  to  I  =  0(  1 ) c | x_  -  x  +  |a2  , 

(ii)  when  c/2  <  c(x)  <  c  for  x  e  (x_-  x+) 

dc  I  x_  -  x+|  <  0(10^,10.)  -  afio^io^.)  <  2dc|x_  -  x+|  , 

when  c/2  <  -c(x)  <  c  for  x  e  (x_,x+) 

dc  |  x_  -  x+|  <  0(10^,00^)  -  afto^.io.)  <  2dc|x_  -  x+| 

Proof.  We  will  prove  the  lemma  when  o  >  0  and  c/2  <  c(x)  <  c,  x_  <  x  <  x+i  other  cases 

are  proved  similarly.  Let  io^x)  and  ui^tx),  x_  <  x  <  x+,  respectively,  be  the  solution 
of  (1.4)  with  io. (x  )  =  io  ,  10  (x_)  =  10  ,  to  (x  )  =  to  and  oo  (x  )  ■  io  .  We  will  denote  by 

P^(x)  the  pressure  for  the  state  oj^(x)  etc.  Clearly  we  have 
Xj(big(x))  >  <1  >  \i(ior(x)  ) ,  x  <  x  <  x+.  We  have  by  integrating  (1.4) 


D.U,  =  Kp  U,  p  U 

£  t  It  r  r 


(3.2)2 


Piui  +  Pl  =  K(pfUi  +  Pi’  +  C  K(y)c(v)Pi(Y)dY 
_ 2  2  x+ 

pu  +p  =K(pu  +  o  1  +  f  K(y)c(v)p(y)dv 
r  r  r  rrr  x  r 


-o- 


13.2)3 


p.E  u  +  p  u  “  K(p.E  u  +  P.u.)  ' 

*  t  *  t  l  £  t  l  X>  K 


pEu  +  p  u  •  K(  p  E  u  +  p  u  )  , 

r  r  r  r  r  r  r  r  r  r 


K(y)  =  exp[-/^  c(z)dz)  , 


K  S  K(x+) 


We  choose  u>  on  S^u^)  such  that 


(3.3) 


<j(u).,H))  ”  0(u).,U  )  =  o  • 
t  It 


It  follows  from  (3.3)  and  the  jump  condition  for  tui^.ii^)  and  (u>,  ,u) 


(3.4),  cfp^  -  pr>  »  P4 -  prur.  o(p4  -  p)  -  Ptut  -  pu  , 


(3.4). 


0(pfuf  -  prV  ’  ( °tue  *  V  '  (prur  +  V  ' 


3<PtUi  ”  °U*  “  <PlUl  +  P£*  '  {P“2  +  P>  ' 


(3.4). 


0(pfEf  '  prV  "  <peVi  +  pfV  ’  tprVr  +  prV 

°(PiEi  *  PE*  *  *°£E£Ui  +  P£U£*  “  * pEu  +  pu) 


We  have  from  (3.2),  and  the  first  identity  in  (3.4),,  that 


*‘pt  ‘  pr’  ’  plut  "  prUr 


Thi3,  and  the  second  identity  in  (3.4>2  yield 


(3.5) 


prUr  *  pu  ”  P<<P,  ~  P>  ”  K(Pj  -  Pr>l 


that 
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I 


1 


we  have  from  (3.2)3  and  the  first  identity  in  0.4)3  that 


_ 2  _ 2  x+ 

(prUr  +  pr}  "  <pi“i  +  ’V  “  0<prUr  “  piU4)  +  4  K<y>c(y)  <pr<y>  "  P^yH^y 


This,  (3.2)1f  and  the  second  identity  in  ( 3.4 ) 3  yield 


(3.5)2  <Pr«£  +  Pr>  “  (P“2  +  P>  -  o(Prur  "  Pu)  +  /x+  K(y)c(y)  (pr<y)  -  p^fyHdy 


From  (3.4)  .j,  we  have 


Pj( o  -  u^)  -  pr(a  -  uf)  =  A  ,  Pl(o  -  u^)  -  p(o  -  u)  i  B 


whence  we  have  from  (3.4) 3  that 


(3.2)1  and  (3.2)3  imply 


E*  ’  Er  +  <piui  '  PrUr)/R  ' 


Ei  ’  E  +  <PiUi  “  pu)/B 


E  +  p/p  ”  E  +  p/p 
r  r  r  r  r  r 


Et +  p/pt  ■ E* +  p/pi 


Thus  we  have  from  the  above  four  identities  that 


(E  +  p  /p  )  -  (E  -  p/p) 
r  r  r 


!  Er  +  Pr/pr  “  Ej  +  (pfcUi  ”  pu,/B  +  p/o 
Ei  ’  (piui  -  Prur)/A  +  pr/pr  _  E*  ‘  p4/pi  +  pt/p4 


+  <P4U£  ~  P“/B  +  P/P 


which  can  be  simplified  to 


] 

i 


( 3  *  S )  3 


p— P  P  — p 

(i,  ♦  ?/».’  -  <<  -  ;/;>  -  t-c-*-  -  zfcrs1 

(JjjtUjj-O)  Kt  l 


He  now  finish  the  proof  of  the  lemma  based  on  (3.5)  ,,  (3.5)2  and  (3.5)3.  He  have  from 


(3.5)1  that 


K  -  1  +  0<1)c|x+  -  x_| 

prur  ”  pu  =  ottPj  “  Pt)  +  <Pr  "  P>)  +  od  -  K)(p1  -  p2)  . 

On  the  other  hand  we  have  from  (3.1)  that 

lpt  "  P4!  “  0(1)c(X4)_,|x+  -  x_| 

lpr  -  pl  -  0(1)c(X(>)'1|x+  -  x_| 

where  X#  is  the  minimum  of  X^tu)  for  all  u  under  consideration.  Since  o  and 
|c(x)|  are  small  compared  to  the  strength  u  of  the  shock  wave  (m^.ti^)  we  know  that 
is  of  the  same  order  of  a  and  it  follows  from  the  last  four  estimates  that 


pfur  -  pu  -  0 ( 1 ) oc | x+  -  x_|(a  +  a) 


(3.6), 


0(1)aclx+  -  x_|a  . 


Note  also  that  pr(y>  -  p^(y),  x_  <  y  <  x+,  is  of  the  same  order  as  a  and  so  (3.6),  and 
( 3.5 ) j  yield 


2ca|x_  -  x^l  ♦  0(1)o2c|x+  -  x_|  a  1  ? 


(pfu2  +  p^)  -  ( pu2  +  p)  >  -j  ca|x_  -  x+|  +  0(1)o2c|x+  -  x_|a  1 

where  we  have  defined  the  strenqth  a  by  the  difference  in  o  just  for  convenience. 
Since  a  is  small  compared  to  a  the  above  estimate  becomes 

3co|x_  -  x+!  >  ( pru 2  ♦  p  )  -  ( pu2  +  p) 

( 3.6 ) 2  , 

>  -  ca|x_  -  x+|  . 

Similarly  we  have  from  ( 3.5)3  that 


(3.6>3 


(Ef  +  —)  -  (E  ♦  £)  -  0(1)co|x_  -  x+|a  1 


Let  us  denote  by  r  any  curve  which  is  the  interaction  of  surfaces  pu“constant  and 
E  +  p/p=*constant.  The  tangent  Tj,  to  r  is  the  cross  product  of 
n^  i  Vpu  and  n2  =  V(E  +  p/p).  Direct  calculations  yield 


(3.7)  T  -  +  X^O,  -1,  -u  -  X^l  +  pev)) 

—  2 

(3.8)  nl  *  n2  "  ”*>eUV  **  ® 

2  S  ^3 

(3.9)  T  •  V(pu  +  p)  - - . 

1+p  v 
e 

Estimates  (3.6).,,  (3.6)2  and  (3.8)  imply  that  there  exists  u  on  r(u>)  Buch  that 

(3.10)  lu  -  o^l  »  0(1)co|x_  -  x+|a  1  . 

From  our  hypothesis  we  know  that  aa  1  is  small  compared  to  a  and  I  X ^  ( <u)  I  is  of  the 
same  order  as  o  for  all  u  under  consideration.  Thus  we  have  from  (3.10)  that 


+  0(  1 )  co|x  -  x+|a 


and  so  (3.9)  yields 

•  y  — .  —2  —  — 2  —  —  —  —  1 

d^a  I  ( Prur  ♦  P  )  _  ( pu  +  p)  |  <  lui  -  u)l  +  O(1)co|x_  -  x+|a  < 

_  y  _  _ 2  —  — 2  — 

2d,  a  I  ( pryr  +  pr>  -  ( pu  +  p)  I 

for  some  positive  constant  d,  depending  only  on  (1.2).  The  above  estimate  and  (3.6>2 
imply  that  for  some  positive  constant  d„ 
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< 3 •  1 1  > ^  0  <  dQ  <  lu  -  ul  (c|x_  -x+|  )  1  <  2dp  <  ••  . 

It  is  clear  from  (2.6)  and  (3.7)  that 

Tj.  •  ^»  >  0 

and  so  (3.9),  (3.6)2  and  (3.11)  imply  that 

(3. 11 ) j  p^  -  p  >  0  . 

It  follows  from  (3.7),  (3.10)  and  (3.11)  that  there  exists  u  on  r^oi)  such  that 

(3.12) ^  lu)  -  o)rl  “  lu>  -  ula  +  Id)  -  u^l 

-  0(l)c|x_x+| {aa  1  +  a2) 

(3. 12 )  2  0  <  dQ  <  lu  -  o)l(c|x_-  x+|  1  <  2dp  <  » 

( 3. 12) 3  p  -  p  >  0  . 

It  is  well  known  that  ( ujg )  is  tangent  to  ( Ug )  up  to  second  order  (c.f.  [2]).  Since 
0)  is  on  S^lo)^),  there  exists  u„  on  S^(u^)  such  that 

(3.13)  lu.  -  ul  ■  0(1)<j2lu  -  ul 


which  and  ( 3. 1 2 ) ^  yield 

lu>r  -  0),  I  4  loir  -  ul  +  lu  -  o>.  I 

(3.14)  «  0(1)c|x_  -  x+|  {era  1  +  a1  +  a2} 

"  0(1)x|x_  -  x+|a2 

This  proves  (1)  of  the  lemmaj  (ii)  follows  from  (3.12)2  and  ( 3. 12 ) 3> 


Q*  E  •  0  e 
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§4.  Numerical  schemes. 

Since  we  will  be  dealing  with  instability  phenomena,  it  is  important  that  numerical 
schemes  used  to  calculate  the  solution  are  effective  so  that  approximate  solutions  do  not 
oscillate  around  the  unstable  solution  or  bifurcate  into  wrong  branch  of  stable  solutions. 
He  will  describe  two  variations  of  the  random  choice  scheme  which  was  introduced  in  [8] , 
simplified  in  [9]  and  is  a  generalization  of  the  Glimm  scheme  for  conservation  laws,  [3], 
[7].  For  convenience  we  write  (1.1)  as 


(1.1) 


,  3f  ( co) 
3t  3x 


g(x,u>) 


and  the  initial  data,  Euler  equations  and  steady  state  equations  become  u(x,0) 


“o(x)' 


(1.2) 


3(i>  3f  (  hi) 
3t  3x 


(1.4) 


3f  ( in) 
3x 


g(x,w)  . 


He  now  describe  the  first  scheme  for  calculating  the  solution  for  the  initial  value  problem 

for  (1.1).  Choose  an  equidistributed  sequence  {a.  }™  in  (0,1),  and  mesh  length 

K  K*0 

Ax  ■  r.  At  -  s  satisfying  the  usual  courant-Friedrichs-Lewy  condition.  The  approximate 
solution  u^fx.t)  ;  uMx.t:  {a^})  is  defined  inductively  as  follows:  Suppose  that 
o)r(x,t)  has  been  defined  for  0  <  t  <  ks,  k  nonnegative  integer.  Then  ur(x,ks+0)  is 
defined  according  to  a^: 

u>f(hr  +  \r‘  hs  +  0)  =  w^Oir  +  a^r,  h  integer  , 

ui^tx,  ks  +  0)  is  a  solution  for  (1.4),  (h-1)r  <  x  <  hr,  h  interger  . 

Thus  ur(x,  ks  +  0)  consists  of  standing  waves  with  possile  discontinuities  at  x  •  hr, 
h  integer.  He  now  resolve  these  discontinuities  so  that  ur<x,t)  is  defined  for 
ks  <  0  <  (k+1)s.  To  resolve  the  discontinuity,  for  example,  at  x  »  0,  we  first  solve 
the  kiemann  problem  for  (1.2)  with 


(ur(x-0,  ks  +  0)  x  <  0 

w(x, 0)  » 

<  01^  ( x+0 ,  ks  +  0 )  x  <  0 

and  denote  the  solution  by  u),(x,t)  «  iji(x/t).  Then  we  perturb  i|i{x/t)  by  (1.4)  to  obtain 
an  approximate  solution  to  (1.1): 

u>(x,t)  =  ^i(£)  when  £  »  x/t  and  t  ♦  0,  and 

bi(x,x/£)  satisfies  (1.4). 

Note  that  <ji(x,t)  exists  locally  in  time  and  is  uniquely  determined  by  u^(x,ks)  provided 
that  (1.4)  can  be  solved  locally,  which  is  the  case  if  X^(u),  i  -  1,3,  is  nonzero  for  all 
states  u  under  consideration.  Finally,  near  (0,ka)  uMx,t)  is  defined  as  a  translation 
of  u(x,t):  ur(x,t  -  ks)  -  w(x,t).  The  discontinuities  (hr,ks)  in  uMx.ks)  are 
resolved  in  a  similar  way.  This  defines  the  approximate  solution  u^Ix.t)  for  all 
(x,t).  Note  that  in  addition  to  solving  the  Riemann  problem,  which  is  needed  in  Glimm 
scheme  for  conservatin  laws,  we  only  have  to  solve  the  ordinary  differential  equation  (1.4) 
at  most  once  for  each  time  and  space  step.  When  a^r/s  is  between  the  maximum  (minimum) 

speed  of  waves  issued  from  (hr,(k-1)s),  ( ( (h+1 )r,(k-1  )s) )  then  ur(x,(k+1)s  +  0)  » 

b>r(x,(k+1  )s  +  0)  for  hr  <  x  <  (h+1)r  and  we  don't  even  have  to  solve  (1.4).  When 
a^r/s  is,  say,  less  than  the  maximum  speed  of  the  waves  issued  from  (hr,ks),  then 
ur(hr+0,(k+1)s  +  0)  is  the  same  as  i|i(a^r/s)  where  is  the  solution  of  the 
corresponding  Riemann  problem  at  (hr,ks).  However  one  has  to  solve  (1.4)  with 

uMhr+0,  ( k+1 )  s  +  0)  as  given  to  find  the  value  of  (h+1  )r-0,  (k+1  )s) . 

The  above  scheme,  when  system  (1.1)  is  reduced  to  conservation  laws  (1.3),  does  not 
zigzag  the  grid  points  (hr.ks)  as  was  done  originally  in  Glimm  13).  This  makes  the 
scheme  numerically  more  efficient.  For  the  present  study  it  is  particularly  convenient, 
since  waves  with  positive  (negative)  speed  are  never  moved  to  the  left  (right).  As  a 
consequence  the  scheme  would  approximate  the  asymptotic  state  of  an  unstable  solution 
accurately,  (see  next  Section).  When  one  is  interested  in  flows  with  a  single  relatively 
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strong  shock  wave  the  following  variation  has  the  advantage  in  tracing  the  evolution  of  the 
shock  wave.  Suppose  that  at  time  t  ”  ks,  the  position  of  the  shock  wave  is  x  ■  x^ , 
(h-1)r  <  xjc  <  hr.  Then  we  disregard  the  grid  points  ((h-1)r,ks)  and  (hr,ks)  and  extend 
the  standing  wave  in  u^tx.ks  +  0)  for  hr  <  x  <  (h+1)r  (or  (h-2)r  <  x  <  (h-1)r)  to 
interval  xk  <  x  <  hr  (or  (h-1)r  <  x  <  Xy )  so  that  ur(x,ks  +  0)  is  discontinuous  at 
x  »  xk  and  continuous  at  x  -  (h-1)r  and  x  -  hr.  This  has  the  advantage  that  the 
position  of  the  shock  wave  is  a  continuous  function  of  time.  In  the  present  study  we  will 
use  the  first  scheme  just  for  simplicity  in  the  presentation. 

The  existence  theory  is  based  on  the  estimates  on  the  total  variation  of  the 
approximation  solution.  This  is  usually  proved  by  induction.  For  later  purpose  we  now 
introduce  some  of  the  notations.  An  I-curve  J  is  a  spacelike  curve  consisting  of  lines 
between  (h+ak)r,ks)  and  (hr,  (k+1)s)  or  between  ((h+a^Kks)  and  ((h+1)r,  (kil)s), 
h  any  integer,  k  any  positive  integer.  J2  is  called  an  immediate  successor  of 
if  J2  lies  toward  a  larger  time  than  J ,  and  they  pass  through  some  grid  points  except 
one.  The  region  between  J1  and  J2  is  called  a  diamond  A.  The  strength  of  waves 
issued  from  (hr,ks)  is  measured  by  its  initial  strength  at  (hr,ks).  Suppose  A  is 
centered  at  ( hr , (k+1 )s) .  Then  waves  entering  A,  besides  those  issued  from  (hr,ks), 
come  either  from  ((h-1)r,ks)  or  ((h+1)r,ks).  The  strength  of  waves  issued  from 
(hr,ks)  is  unchanged  at  (hr , (k+1 )s-0 ) .  Those  issued  from  (h-1)r,ks)  (or  (h+1)r,ks)) 
have  positive  (or  negative)  speed  and  the  difference  of  their  strength  between  times 
ks  +  0  and  (k+1)s-0  is  due  to  the  geometry  of  the  duct  between  x  »  (h-1)r  and  x  «  hr 
(or  between  x  «  hr  and  x  -  (h+1)r)  and  is  estimated  by  Lemma  3.1.  The  waves  leaving 
A  are  then  the  result  of  nonlinear  interaction  for  (1.2)  of  waves  inside  A  at  time 
(k+1)s-0  just  mentioned.  Such  interactions  have  been  studied  in  [3)  and  [6). 
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§5.  Unstable  flows. 

In  this  section  we  study  transonic  flows  with  positive  speed  along  a  contracting 
duct.  For  definiteness  we  assumet 


(5.1) 


c( x)  =  -  >0  for  0  <  x  <  1 

a(x) 

c(x)  *0  otherwise  . 


Consider  the  solution  ai„(x,t  >  -  « *(x)  of  (1.1)  which  consists  of  standing  waves  and  a 
standing  shock  wave  at  x  «  x» ,  0  <  xA  <  1 , : 


(5.2) 


(x) 

0  «  X  <  X, 

Wjfx) 

X#  <  X  < 

(!>(-«•) 

x  <  0 

lil(+*) 

X  >  1 

where  u^lx)  and  <»2(x)  are  solutions  of  (1.4),  (u^  (xt)  ,a>2<xt) )  is  a  1->shock  wave  with 


(5.3) 


o(tal1  (x,)  ,<1L  (x#)  )  -  0 


and  tii(-<»)  =  <^(0),  o»(+«)  S  u2(1).  Before  we  study  the  instability  of  standing  shock 
waves,  we  fist  show  that  wave  patterns  of  other  kinds  are  stable.  We  will  then  show  that  a 
small  perturbation  of  a  standing  shock  wave  causes  the  unstable  flow  to  approach  one  of  the 
stable  flows.  By  a  nonlnteractlng  wave  pattern  we  mean  a  solution  of  (1.1)  consisting  of 
elementary  waves  for  Euler  equations  (1.3)  in  the  region  x  t  (0,1)  and  of  standing  waves 
and  a  possible  standing  shock  wave  in  x  e  (0,1).  Moreover,  these  waves  do  not  interact  in 
the  positive  time  direction,  in  other  words,  they  do  not  occupy  the  same  space  position  for 
any  positive  time. 

Theorem  5.1.  Suppose  that  u(x,t)  is  a  noninteracting  wave  pattern  for  (1.1)  with  the 
property  that  u(x,t)  does  not  contain  any  standing  shock  wave  and  all  waves  In  u(x,t) 
are  weak  and  X^u)  j*  0  for  u  «  u(x,t).  Then  w(x,t)  is  asymptotically  stable  in  the 
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sense  that  if  an  initial  data  <o(  x,0 )  is  such  that  the  total  variation  of 
to(x,0)  -  to(x,0)  is  small,  then  the  solution  w(x,t)  of  (1.1)  tends  to  an  asymptotic 
state  consisting  of  elementary  waves  which  are  close  to  those  in  to(x,t)  both  in  position 
and  strength. 

Proof:  We  may  let  the  perturbation  u(x,0)  -  to(x, 0)  be  so  small  that  u(x,t)  is  never 
sonic  and  so  the  schemes  described  in  the  previous  section  may  be  applied.  Note  that  since 
<n(x,t)  contains  no  standing  shock  wave,  the  standing  wave  in  it  is  either  strictly 
supersonic  or  strictly  subsonic.  Although  Liu  [8]  studies  only  flows  which  are  supersonic 
(or  subsonic)  throughout  the  duct,  the  ideas  employed  there  can  be  generalized  to  prove  the 
present  theorem.  We  omit  the  details. 

Q.  E.  D . 

Theorem  5.2.  The  wave  pattern  ui#(x,t),  (5.2)  and  (5.3),  is  dynamically  unstable  provided 

that  |c(x)|  is  small  as  compared  to  the  strength  aQ  of  the  standing  wave 
( <u  (x  ) ,  to  (x  ) ) .  More  precisely,  there  exists  (to(x,0)  such  that  the  total  variation  of 

i  *  2  * 

bit x , 0 )  -  to,  ( x ,  0 )  is  arbitrarily  small  and  as  time  goes  to  infinity  co(x,t)  tends  to  one 
of  the  stable  asymptotic  states  without  a  standing  shock  wave  as  studied  in  the  previous 
theorem. 

Proof:  For  definiteness  we  assume  that  0  <  x,  <  1.  We  choose  the  perturbation  to(x,0) 
to  consist  of  standing  waves  and  a  1 -shock  wave  at  x  **  x,  with  positive  speed.  This  can 
be  done  quite  easily  by  choosing  a  state  to  on  Sj(<0j(x,))  which  is  between  to^tx,)  and 
to^lx,)  and  arbitrarily  close  to  u^(x,  ),  and  define  to(x,0)  by: 

to^tx),  0  <  x  <  x, 

<o(x,0)  =  i^(x)  x,  <  x  <  1 

u^(x,t)  otherewise 

where  uij(x)  is  a  standing  wave  with  <^(x,)  ■  to.  u>(x,0)  is  a  small  perturbation  of 
<o,(x,t)  since  to  is  close  to  ^(x,).  Moreover  we  have 


o(u(x#  -  0,  0),  u(x#  +  0,  0))  >  0  . 

For  convenience  we  choose  Ax  -  r  such  that  ( x* ,  0 )  is  a  grid  point.  He  will  use  the 

first  scheme  described  in  the  previous  section.  He  first  establish  the  existence  of  the 

solution  u(x,t).  The  first  step  is  to  establish  a  uniform  bound  on  the  total  variation  of 

the  approximate  solution  Uj(x,t) .  He  will  show  inductively  that  there  is  a  relatively 

strong  1-shock  wave  in  ur(x,t)  which  originates  at  (x*,0)  and  persists  for  all  time. 

Suppose  that  it  crosses  an  I-curve  J  at  x  «  hQr  and  has  strength  a* .  Then  we  define  a 

nonlinear  functional  F(J)  as  follows: 

F(J)  =  L(J)  +  K{Q,(J)  +  Q(J)  +  Q  (J)  +  Q.  (J)  }  , 

1  2  <5  Q 

L( J )  5  J { |  oil  :  a  the  strength  of  any  elementary  wave  crossing  J)  , 

Q,  (J)  i  qJ(J)  +  Q*(J)  +  , 

Q*(J)  =  J { 1  oil  /*  |c(x)|dx  :  a  the  strength  of  any  1-wave  crossing  J  and  issued 

from  (hr,ks)  with  h  <  hQ}  , 

C>![(J)  S  T{lal  f"  |c(x)  I  dx  +  I  al  f?’r  |c(x)|dx:  a  the  strength  of  any  1-wave 
1  n^r  a  r 

0  o 

crossing  J  and  issued  from  (hr,ks)  with  h  >  hQ}  , 

5  £  r  !<=(>*)  lax  , 

q®(J)  s  £{|oBl  :  a  and  B  strengths  of  1-waves  crossing  J  and  at  least  one  of 

them  is  a  shock  wave)  , 

Qj(J)  £  y  { |  ci I  fhr  :  a  the  strength  of  any  2-wave  crossing  J  and  Issued  from 

(hr.ksl)  , 

Q3(J)  H  Q^(J)  +  Q*(J)  . 

Q'(J)  5  ^{lal  f"  |c(x)|dx  :  a  the  strength  of  any  3-wave  crossing  J  and  issued 
3  ■  hr 

from  (hr,ks>)  , 

Q® ( J )  =  y { 1 aBl  :  a  and  B  strengths  of  any  two  3-waves  crossing  J  and  at  least 

one  of  them  is  a  shock  wave)  , 

Qd(j)  ^  y { I a6l  :  a  and  8,  respectively,  strengths  of  a  i-wave  to  the  right  and  a 

j-wave  to  the  left  crossing  J  and  i  <  j)  , 
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Similarly,  when  h  -  h2,  waves  in  have  positive  speed  and  more  from  x  -  (h-1)r 

to  x  •  hr  during  the  time  t  e  (ks,(k+1)s).  We  set 

QU>  H  (^h-1)r|c(x,|dx)  *  |  '‘VVi1  +  l(“l'"!2,i<“s'“4,jl 

+  J{l(w1.<n2)1<(a3,w4)il  :  at  least  one  of 
and  (uij/Ug)^  is  a  shock  wave}. 

With  this  definition  of  Q(A)  we  have  from  Lemmas  3.1  and  3.2  and  results  on  the 
interaction  of  elementary  waves  for  conservation  laws,  [3],  [8],  that 


a.  *  +  (u3,a>4)l  *  Ofl)C(A) 


where  is  the  strength  of  the  i-wave  leaving  A. 

Suppose  that  (S.4)  and  (5.5)  have  been  shown  for  an  J  «  ](  and  J2  is  an  immediate 

successor  of  Jj,  A  the  diamond  between  them.  According  to  <5.6),  waves  may  change  their 

strength  only  due  to  linear  combining  or  cancelling,  and  nonlinear  interation  as  measured 

by  Q.  Since  the  only  wave  crossing  0  is  the  1 -shock  wave  with  strength  a  and  issued 

from  (x*,0),  it  follows  that  the  strength  of  this  shock  wave  as  it  crosses  J,  is  of  the 
strength  +  o(1)Q(  and  all  other  waves  have  a  total  strength  0(1 )Q(  1),  where 
Q(  1 )  is  the  region  below  Jr  It  follows  from  the  induction  hypothesis  (5.5)  that 


*1 

<  2  (Q CO)  -  Q(J)1  <  2  fx'  |c(x)  |dx  , 


where  x  ■  x,  is  the  position  of  the  relative  shock  wave  as  it  crosses  J^.  Thus,  in 
particular,  the  relative  strong  shock  wave  has  strength 


xl  “n 

a.  -  a  (1  +  0(1)  /  1  |c(x)  |dxl  >  t2- 


as  it  crosses  J..  The  speed  of  this  shock  wave,  according  to  Lemma  3.2  and  (5.7),  is 


(5.9) 


>  o(x.)  +  V2  d  I c ( y )  I <Jy  +  0(1  )Q ( A1 ) 

1  X1 

>  o(xt)  ■*•  -  d  /  | c { y )  | dy  . 

j  x# 

Note  that  o(x#)  la  assumed  to  be  positive  (though  arbitrarily  small)  and  a  wave  with 
positive  speed  may  travel  only  to  the  right  according  to  the  scheme  introduced  in  the 
previous  section.  Consequently.  x4  >  x^  and  (5.9)  implies  by  induction  that  the  relative 
strong  shock  wave  always  have  positive  speed  at  least  up  to  the  point  when  it  crosses 
we  now  establish  (5.4)  and  (5.5)  for  J  “  Jj.  Baaed  on  the  above  description  of  the 
relative  strong  1-shock  wave  we  know  that  1-waves  to  the  right  (left)  of  this  shock  wave 
have  negative  (positive)  speed.  This,  along  with  estimate  (5.6),  imply  that 

F(J2)  -  FfJ^  <  -  V2Q(A)  <  0 
Q(J2>  -  Q(Jt>  <  -  V2QU) 

whence  we  have  easily  (5.4)  and  (5.5)  for  J  “  J2>  Actually,  the  functionals  F(J)  and 
Q(J)  were  defined  so  that  these  two  crucial  estimates  hold  (see  also  (BJ).  Details  are 
left  to  the  reader. 

To  finish  the  proof  of  the  theorem,  we  show  that  the  relatively  strong  shock  wave 
leave  the  region  x  e  (0,1)  in  finite  time  so  that  the  solution  becomes  a  small 
perturbation  of  a  stable  wave  pattern  as  studied  in  the  previous  theorem.  We  know  from 
(5.9)  that  the  speed  of  the  shock  wave  is  always  positive.  It  remains  to  show  that  as  a 
result  of  the  equidistributedness  of  the  sequence  {a^}  and  (5.9),  the  time  limit  for  the 
shock  to  leave  x  c  (0,1)  is  independent  of  the  mesh  length  Ax  «  r.  Assume  for  the 
moment  that  the  shock  wave  propagate  along  a  continuous  curve  (which  would  be  the  case  if 
the  second  scheme  of  the  previous  section  was  used).  Let  x  »  x(t)  be  the  position  of  the 
shock  at  time  t,  then  x' (t)  “  o(x)  and  we  have  from  (5.9): 
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x“(t)  >  J  d  C ( X ( t ) )  x  * (t)  . 

For  the  approximate  solution,  the  above  estimate  holds  module  an  error  which  depends  on  how 
well  equidistributed  the  sequence  {a^}  is  and  goes  to  zero  as  ix  =  r  tends  to  zero, 
([7]).  Thus  for  simplicity  we  may  argue  with  the  above  inequality.  Since  d  is  a 
Dositive  constant,  c(x)  is  positive  for  0  <  x  <  1,  x(0)  e  (0,1),  x'(0)  is  positive 
(though  it  may  be  arbitrarily  small),  the  inequality  implies  that 

x(t)  >  x'(0)  exp  [j  jjj  c(x(s))ds)dz 

and  so  there  exists  a  finite  time  t«  independent  of  Ax  <■  r  such  that  x(t)  >1  for 
t  >  t# •  In  other  words,  the  relative  strong  shock  wave  leaves  the  region  0  <  x  <  1 
before  finite  time  t».  This  completes  the  proof  of  the  theorem. 

Q.E.D. 


-24- 


$6.  Stable  flows. 

We  already  know  that  a  noninteracting  wave  pattern  containing  no  standing  shock  wave 
Is  stable  and  a  flow  along  a  contracting  duct  which  contains  a  standing  shock  wave  Is 
unstable.  To  complete  the  present  study  It  remains  to  show  that  a  flow  with  a  standing 
shock  wave  along  an  expanding  duct  la  stable.  Since  for  an  expanding  duct,  an  asymptotic 
(l.e.  noninteracting)  wave  pattern  depends  uniquely  and  smoothly  on  Its  end  states  at 
x  ”  ±»,  [101,  It  suffices  to  show  that  a  small  perturbation  of  a  flow  with  a  standing 

shock  wave  tends  to  an  asymptotic  wave  pattern  as  time  tends  to  Infinity.  Throughout  this 
section  we  assume: 


(6.1) 


c(x)  =  -  *  <0  for  0  <  x  <  1  , 

c(x)  >0  otherwise  . 


Theorem  6.1.  Suppose  that  u(x,0)  contains  a  1-shock  wave  (u>_,u+)  at  x#  e  (0,1)  with 

X 1  ( )  <  0  <  Xj(w+),  and  both  |c(x)|  and  |o(u)_,u+)l  are  small  compared  to  the 

strength  a  of  the  wave  (u>_,u+).  Suppose  also  that  w(x,0)  is  a  small  perturbation  in 
total  variation  of  standing  waves  for  -»  <  x  <  x#  and  for  x#  <  x  <  •».  Then  there  exists 
a  solution  u(x,t)  of  (1.1)  which  tends  to  a  noninteracting  wave  pattern  as  t  tends  to 
infinity. 

Proof:  The  difference  between  the  situation  here  and  that  in  Theorem  5.2  is  that  although 
there  still  exists  a  relatively  strong  shock  wave  for  all  time,  this  shock  wave,  unlike  the 
corresponding  one  in  Theorem  5.2,  stays  nearby  the  unperturbed  standing  shock  wave. 
Consequently,  we  have  to  show  that  the  shock  wave  decelerates  in  general  and  so  would  not 
cause  too  strong  nonlinear  interactions  if  it  were  to  oscillate  along  the  expansion  portion 
of  the  duct.  For  this  we  define  a  different  nonlinear  functional  F(J)  as  follows: 

Suppose  that  the  shock  wave  has  speed  a  "  c(t),  strength  a  and  position  x  »  hgr  “ 

hQ(t)r  when  it  crosses  J  (at  time  t).  We  set 

F(J)  =  L(J)  +  K ( J)  +  Q2(J)  +  Q3(J)  +  Sd<J)) 

c 

where  L  and  Q's  are  the  same  as  those  in  the  previous  section  except  that  Q^(J)  is 
set  to  be  zero  and  L(J)  is  redefined  as 
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L( J )  5  Kiel  +  J  lal  i  a  the  strength  of  any  wave  crossing  J}  . 

The  proof  of  uniform  boundedness  of  the  total  variation  in  x  of  the  approximate  solution 
uir(x,t)  is  similar  to  the  proof  of  Theorem  S.2.  The  notable  difference  here  is  that 
since  we  do  not  know  a  priori  the  sign  of  the  speed  of  the  relative  strong  shock  wave,  we 
cannot  foresee  the  potential  interaction  of  the  shock  with  the  geometry  of  the  duct  (which 

(J 

was  taken  care  of  by  in  Theorem  S.2).  Instead,  we  notice  from  Lemma  3.2  that  when  the 

shock  moves  between  x  »  0  and  x  -  1,  where  the  duct  is  not  uniform,  it  tends  to 
decelerate  and  so  |o|  tends  to  decrease.  This  shows  that  the  first  part  of  L(J)  tends 
to  decrease.  This  is  sufficient  to  dominate  ne  increase  in  the  remaining  parts  of  F(J) 
as  the  result  of  the  nonlinear  interaction  of  the  relative  strong  shock  wave  with  the 
duct.  Therefore,  in  place  of  (5.5),  we  have 

(6.2)  Q(A)  <  O(1)o(0)  +  2{Q(0)  -  Q(J))  . 

The  above  statements  are  shown  using  Lemmas  3.1,  3.2  and  analogous  arguments  used  in  the 
proof  of  Theorem  5.2j  details  are  omitted. 

We  now  study  the  asymptotic  behavior  of  the  solution.  Given  any  e  >  0,  it  follows 
from  (6.2)  that  there  exists  T  <*  T(  e)  such  that 

(6.3)  Q{t  >  T}  <  E  . 

By  choosing  M  »  M( e)  sufficiently  large  we  have 

(6.4)  total  variation  x  (u(x,T)  :  |x|  >  M)  <  e  . 

Through  (H,T)  and  ( -M,T)  we  draw  generalized  i-characterlstic  curves  x^  and  x^» 

1  »  1,2,3  which  travel  either  with  i-charaeteristic  speed  or  i-discontinuity  speed  (4), 
[7],  The  region  between  x^  *nd  is  denoted  by  0^ .  Due  to  the  strict  hyperbolicity 

of  the  system,  x*  do  not  intersect  Xj'  i  J*  j ,  after  finite  time  Tg  >  T.  After  time 
Tq,  the  i-waves  outside  are  either  produced  by  interaction  or  issued  from  (x,T), 

I x t  >  M«  and  so  (6.3)  and  (6.4)  imply  that  the  total  amount  of  i-waves  after  time  Tg  and 
outside  is  0(1 )e.  Thus  in  0^ ,  i-waves  cancel  and  combine  and  behave  (mod  O(I)e) 

like  waves  for  scalar  conservation  law,  (Section  10,  [8)).  Consequently,  module  the  error 
0(1 )e  3-waves  tend  either  to  a  3-shock  wave  or  a  3-rarefaction  wave,  2-waves  tend  to  a 
traveling  wave  and  1 -waves  tend  to  a  shock  wave  due  to  the  presence  of  the  relative  strong 
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shock  wave  in  the  solution.  Moreover,  outside  the  1-shock  wave  and  in  the  region 

0  <  x  <  1,  the  solution  equals  (mod  0(1 )e)  standing  waves.  To  show  the  closeness  of  such 

elementary  waves  to  a  noninteracting  wave  pattern  it  remains  to  show  that  when  the  1 -shock 

wave  stays  in  the  region  0  <  x  <  1  its  speed  is  close  to  zero.  In  fact  we  will  show  that 

2+|x0( 

in  this  case  its  speed  is  0(1)/e"  after  time  Tg  +  — -  where  x  «  Xg  is  the  position 

of  the  1-shock  wave  at  time  Tg.  For  this  we  first  note  that  the  speed  o(t)  of  the  shock 
wave  at  time  t  is  a  Lipschitz  function  of  t  (mod  0(1)e)  in  the  following  sense:  By 
Lemma  3.2  we  have  for  t_  >  t.j  >  Tg 


lo(t2)|  <  !  o<  1 1 )  +  0(1)Q{t  c  (t^tg)}  , 


-  ?  1 1  3  f  +  0(1)e  . 

_  _  2+(x  ) 

Thus  if  I  o( t )  1  <  /e  for  some  t  e  (tq,  tq  +  — - - )  then  |<j(t)|  <  ft  +  0(1)e  «  2/e  for 

lx  |  +  1 

t  >  t„.  It  remains  to  show  that  if  lo(t)|  >  ft  for  all  t  e  (t_,  T_  +  - )  then  the 

0  '00  c  ; 

u0|+2 

1 -shock  wave  does  not  stay  in  the  region  0  <  x  <  1  after  time  TQ  +  - — - — .  In  fact  itis 
clear  that  when  the  speed  of  the  one  shook  is  larger  than  /I  (or  less  than  -ft)  for  all 
IxqI+1  IXqI+1 

t  e  (Tq,  Tq  ♦  — - - )  then  at  time  Tg  +  — - -  the  position  of  the  shock  wave  is  to  the 

riqht  of  x  «  1  (or  left  of  x  «  0)  and  by  (6.4)  its  speed  after  that  time  is 
ft  -  0(1  )c  >  Vj  ft  (or  -/c  +  0(1)e<-  ^2  ft)  and  always  stays  in  the  region  x  >  1 
(or  x  <  0)  for  any  larger  time.  This  completes  the  arguments. 

When  the  perturbed  data  u(x,0)  equal  to  constant  states  both  to  the  right  and  to  the 
left  of  a  finite  Interval,  the  convergence  of  the  solution  u(x,t)  to  noninteracting  wave 
patterns  is  of  algebraic  rates.  For  this  more  precise  estimates  are  needed,  we  will  not 
pursue  that  here  (c.f.  (6),  18]).  This  completes  the  proof  of  the  theorem. 
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